Universal routes to spontaneous PT-symmetry breaking 
in non-hermitian quantum systems 



o 

(N 

> 
O 

in 



^ ■ 

Oh 
I 



> 

00 




NxT 



(b) orthogonal 
5A 





unitary 



Henning Schomerus 

Department of Physics, Lancaster University, Lancaster, LAI 4YB, United Kingdom 

(Dated: November 2010) 

■PT-symmetric systems can have a real spectrum even when their Hamihonian is non-hermitian, 
but develop a complex spectrum when the degree of non-hermiticity increases. Here we utilize 
random-matrix theory to show that this spontaneous "PT-symmetry breaking can occur via two 
distinct mechanisms, whose predominance is associated to different universality classes. Present 
optical experiments fall into the orthogonal class, where symmetry-induced level crossings render the 
characteristic absorption rate independent of the coupling strength between the symmetry-related 
parts of the system. 
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Non-hermitian quantum systems generally have a com- 
plex energy spectrum, with imaginary parts of the ener- 
gies related to decay or amplification rates. However, 
when loss and gain are in balance the spectrum can still 
be real. One intensely researched route to try and achieve 
such a balance is to couple two identical systems sym- 
metrically and then induce opposite amounts of gain and 
loss into the two parts, as illustrated in Fig. [Ija) (ll-[l3j|. 
The Hamiltonian then possesses a combined parity (V) 
and time-reversal (T) symmetry, and its secular equation 
is real. However, this does not guarantee a real spec- 
trum; as the level of non-hermiticity (loss and gain) is in- 
creased, pairs of complex-conjugate energy levels appear 
[^-Q. This phenomenon of spontaneous T^T-symmetry 
breaking has gained recent prominence because it leads 
to optical effects such as double refraction, solitons and 
non-reciprocal diffraction patterns, which provide mech- 
anisms for the design of unidirectional couplers and left- 
right sensors 0, Q, concepts that are now being real- 
ized experimentally in a variety of optical settings [^. 
Over the past months, these systems were proposed for 
at-threshold lasers jlC| and laser-absorbers [ll], [13] . In 
turn, these developments have instigated a deeper theo- 
retical understanding of the role of the dynamics (such 
as the consequences of Anderson localization and wave 
chaos 12], as well as interactions [13]). In this paper, we 
establish distinct universality classes which directly affect 
the nature of spontaneous T^T-symmetry breaking. 

To do so, we derive random-matrix ensembles where 
loss and gain in the two parts of the system are im- 
plemented by a uniform rate fi, while coupling is es- 
tablished through N channels with transmission prob- 
ability T; the mean level spacing of the decoupled parts 
is A [15[. We find that the mechanism behind sponta- 
neous PT-symmetry breaking depends on whether the 
hermitian limit fj, — is time-reversal symmetric or 
not, amounting to a predominance of symmetry-induced 
level crossings or level repulsion, as illustrated in Fig. 
(lib). This results in different characteristic scales fipT 
of amplification/absorption governing the transition from 




FIG. 1: (Color onhne) (a) Sketch of a nonhermitian 'PT- 
symmetric system, where a region with absorption rate fi 
(and mean level spacing A, left) is coupled symmetrically via 
a tunnel barrier (supporting A'^ channels with transmission 
probability T) to an amplifying region with a matching am- 
plification rate (right). Below this, the scattering description 
of the system, (b) Two routes to spontaneous PT-symmetry 
breaking, depending on whether the hermitian limit jj, — 
is T-symmetric (orthogonal class displaying level crossings, 
left) or not (unitary class displaying avoided crossings, right). 
Shown are real eigenvalues of a random Hamiltonian H [Eq. 
(Q] as function of T for fixed jj. — (left of dashed line), and 
then as a function of /i for fixed T — 1 (right of dashed line). 
Complex-valued levels (formed by level coalescence at /i > 0) 
are not shown. Here fio ~ \/iVA/27r, and we set A'^ — 10. 



an essentially real to an essentially complex spectrum. 
Present optical experiments and theoretical studies con- 
cern systems without magneto-optical effect (the orthog- 
onal symmetry class), in which the hermitian limit is 
time-reversal symmetric. In this case /ipT ~ ■\/]VA/27r = 
fiQ becomes fully independent of the coupling strength 
as soon as T surpasses a parametrically small threshold 
Tc ^ 1/-/V, thereby exhibiting a level of universality that 
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jugation, 'PT'-symmetry results in the relation [IC 
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FIG. 2: (Color online) Gradient plots of the ensemble- 
averaged fraction /(/x, T) of complex energy levels (among all 
levels witiiin a range of energies over which the mean spacing 
A can be assumed constant), for the orthogonal class (left) 
and the unitary class (right). In (a) and (b), ^ is scaled to 
/xq. The bottom panels show same data with fj, scaled to 
/xo = ^J.o / ^/n^T/NT (c) and fi'^ = Vt^o (d). Numerical 
results with N = 50. 



goes beyond what is normally encountered in mesoscopic 
systems. For weak coupling (T < Tc), ^ipr ^ VNTfiQ. 
Adding magneto-optical effects to the system essentially 
changes the nature of the transition. In this case (the uni- 
tary symmetry class), /ipT ^ VTfJ-o hi the full range of 
weak and strong coupling. These findings are illustrated 
in Fig. El 

Random-matrix ensembles. — To derive the appropriate 
random-matrix ensembles we formulate a quantization 
condition based on scattering theory 3, 12|- The 
N X TV-dimensional scattering matrix 
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of the left subsystem can be expressed in terms of an 
M X M-dimensional Hamiltonian H, which is real sym- 
metric (a member of the standard Gaussian orthogonal 
ensemble, GOE) if the hermitian limit = is time- 
reversal symmetric, and complex hermitian (a member 
of the Gaussian unitary ensemble, GUE) if this is not the 
case [ill . We assume M ^ N ^ 1 and denote the mean 
level spacing in the energy range of interest by A. The 
M X N coupling matrix V then fulfills VV^ — diag(wm), 
where N diagonal entries Vm — AM/tt correspond to 
fully transparent channels, and M — N entries = 
describe the closed channels . Adopting a basis where 
the time-reversal operation T is identical to complex con- 
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for the scattering matrix of the right subsystem. The 
tunnel barrier is described by refiection amplitudes r — 
— Vl ^ T and transmission amplitudes t = ^^/T. As 
shown in the lower part of Fig. [IJa), these scattering 
matrices relate amplitudes of left and right propagating 
waves at the two interfaces of the tunnel barrier. The re- 
quirement of consistency of these relations results in the 
quantization condition 



det 



r t 
t r 



Sl 
Sr 



1 



0, 



(3) 



which can be rearranged into an eigenvalue problem 
det {E — H) = with effective Hamiltonian 
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H -ifi T 
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The positive semi-definite coupling matrix F — diag (7m) 
now incorporates the finite transmission probability of 
the barrier; its N non- vanishing entries read 7™ = 
[VT/il + VT^)]AM/Tr = 7 [ij- 

Numerical evaluation. — Before engaging in an analyt- 
ical discussion of the different routes to spontaneous VT- 
symmetry breaking [illustrated in Fig. [ijb)], we put for- 
ward numerical results which illustrate the physical con- 
sequences of the points to be made below. These results, 
presented as (color) gradient plots in Fig. [51 concern the 
fraction /(/x, T) of complex- valued energy levels within a 
range where the mean level spacing can be assumed con- 
stant. We fix M = 1000, which ensures a large number 
of levels within the range in question, and set iV — 50 

Figure [21^ a) shows results for the orthogonal ensemble, 
with /i scaled to (Mq — ^/NA/2n. We see that above a 
small threshold Tc (to be determined below as Tc ^ ^/N), 
the transition from the real spectrum (/ = 0, obtained 
for fi — 0) to a spectrum which is partially real and par- 
tially complex spectrum {f 1/2) indeed occurs on the 
scale fipT ~ fJ-Q, and then is independent of the value of 
T. Only for T < Tc, ^ipT ^ VNTfiQ = fix is coupling- 
dependent. In order to get a unified view over both 
regimes, we plot in panel (c) the same data, but with 
fi scaled to = ^o/\/l + 1/NT, which interpolates be- 
tween fiT for T <^ Tc and /xo for T ^ Tc. The conver- 
gence of gradient lines for T — >■ indicates that for weak 
coupling the transition becomes more abrupt. 

Figure [2{b) shows the corresponding results for the 
unitary ensemble, where fi is again scaled to ^o- Here 
we find that a systematic T-dependence persists across 
the full range of coupling strengths. As shown in panel 
(d), this dependence takes the form /ipT ~ n/T/zq = fi'x- 
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Now, the only difference between the strong and weak 
coupling regimes is a factor of order 1 . In further contrast 
to the orthogonal case, for weak coupling the transition 
remains smooth; however, since /i^ = iit/VN ^ ^j-, it 
then occurs at a far smaller deviation from hermiticity. 

Underlying mechanisms. We now show that the fea- 
tures reported above originate from two distinct mecha- 
nisms of spontaneous T'T-symmetry breaking. 

It is instructive to start in a regime which can be 
treated perturbatively. For n — 0, the effective Hamilto- 
nian H [Eq. is hermitian and all its eigenvalues are 
real. For T = (F = 0), on the other hand, the spectrum 
is a superposition of two level sequences = ± i/i, 
which are all complex if /i ^ 0; here e^. are the eigen- 
values of H. Therefore, in regard to the question of how 
many levels are complex, the limits T, /i — ;> do not com- 
mute. Nonetheless, for T = = the spectrum reduces 
to the superposition of two degenerate level sequences 
e/c, so that quasi-degenerate perturbation theory applies. 

(k) 

Denote by ipm the wave function of H corresponding 
to eigenvalue e^; in random-matrix theory, this is a ran- 
dom normalized vector with average |'(/'m''p = 1/M. Re- 
duced to the symmetric and antisymmetric extension of 
this wavefunction across the whole system, the effective 
Hamiltonian takes the form 
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whose eigenvalues become complex for /ipT 
I Sm[''/'™^]^7™l- Therefore, on average (and using 7 
VTAM/2n for T < 1) 



MPT 



NVTA/2TT 



jiT (orthogonal class), 
pi'rp (unitary class). 



(6) 



which recovers the numerical scales in the weak coupling 
regime. 

Note that the two expressions for ^pt differ by the 
parametrically large factor ^ ^/N. Mathematically, this 
arises because Vm^ is real in the orthogonal class and 
complex in the unitary class; physically, it amounts to 
vastly different tunnel splittings. This difference signifies 
that in the orthogonal class, the levels of the originally 
degenerate sequence £k from the two subsystems quickly 
cross as T is increased. A second route to PT-symmetry 
breaking then becomes available, which involves two en- 
ergy levels that are non-degenerate for T = 0. In order to 
describe this case we reformulate the problem by starting 
with /i = 0, and exploit the thus-emerging P-symmetry 
in the orthogonal class to transform the effective Hamil- 
tonian to 



Hv = 
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We denote by the two level sequences of iJ±F. Since 
F is positive semidefinite these sequences arise from the 



sequence Sk by an oppositive shift which is approximately 
rigid. From the resulting combined sequence, consider 
two levels and which lie adjacent to each other; the 
corresponding eigenvectors are ijj^''^^ and Finite 
fi mixes these levels, which is embodied in the reduced 
Hamiltonian 



in{i^'-''+^^'-^- 



(8) 



and — sequence, {xp^'^^'^ 



Now, treating 2F as a perturbation which connects the -I- 

('-)\ « Because 

je^ — £,^| = ©(A) is small compared to the shift due 
to the coupling, the denominator can be estimated as 
4 - 4 ~ - 4 ~ -(V'('+^|2F|i/'('+)). The coupling 
strength drops out, and on average ^ 
1 /N, i.e., the mixing is small. As a result, the level pair in 
question becomes complex for /z^|(V'^'~^|V'*-'^^-')P ~ i^k~ 
eif ^ A2, i.e.. 



^ VNA/2TT = fiQ (orthogonal class, T > l/N). 

(9) 

As indicated, comparison of this expression with Eq. ([6]) 
implies that this mechanism becomes favorable around 
T = Tc^ l/N. 

This analysis of strong coupling does not apply to the 
unitary case, which does not display 7^-symmetry for /i = 
0. In the T'-basis, in place of Eq. ^ we then have 



Hv = 



RcH + T ilmH + ifi 
ilmH + ifl RcH - F 



(10) 



Consequently, finite coupling not only results in a far 
reduced systematic shift of the levels, but also in a direct 
mixing of levels in the individual sequences. Therefore, 
instead of level crossings one encounters level repulsion. 
This difference is illustrated in Fig.[ljb), which shows the 
evolution of energy levels as T is increased from to 1 
(while /i = 0) , and the subsequent fate of real levels as fi is 
increased from to 4^o (while T = 1); pairwise coalescing 
levels become complex, and then are no longer shown. 
In the orthogonal class, such pairs trace back to well- 
separated levels , from the two different sequences 
(which are distinguished by the opposite slopes of the 
levels for increasing coupling). In contrast, in the unitary 
class the coalescing levels trace back to originally closely 
spaced or degenerate levels, even when the coupling is 
strong. 

Finally, we point to an alternative scenario where the 
coupling- independent fiPT ~ Mo becomes relevant even in 
the unitary symmetry class. Observe that by definition, 
for a hcrmitian Hamiltonian the T operation (complex 
conjugation) is equivalent to transposition, an operation 
that we denote by T' ■ However, for non- hermitian sys- 
tems there is a physical difference: Absorption and am- 
plification, taken by themselves, break T-symmetry, but 
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preserve T'-symmetry; the latter is broken by magneto- 
optical effects. We find that a combined PTT'-symmetry 
still results in a spectrum with levels that are either real 
or occur in complex conjugate pairs. Compared to the 
case of 'PT'-symmetry, experimental implementation sim- 
ply requires to invert the magneto-optical effects in one 
part of the system. The effective Hamiltonian then takes 
the form 

which differs from (j4|) when H is complex (i.e., in the uni- 
tary symmetry class). In the parity basis, the Hamilto- 
nian takes the form of Eq. ([7]) even for unitary symmetry. 
Coupling now induces level crossings, and the transition 
to the complex spectrum is governed by the same charac- 
teristic scales fj,T and as encountered in the orthogonal 
symmetry class of PT-symmetric systems. 

Conclusions. — In summary, we identified two routes to 
the formation of complex energy levels in non-hermitian 
quantum systems with PT-symmetry (spontaneous VT- 
symmetry breaking). The predominant mechanism de- 
pends on whether or not the hermitian limit possesses 
time-reversal symmetry (orthogonal or unitary universal- 
ity class, respectively). Present optical experiments fall 
into the orthogonal class, where level crossings result in 
a characteristic absorption/amplification rate ^J-'PT which 
is independent of the coupling between the symmetry- 
related parts of the system (unless the coupling is very 
weak). The unitary class features strong level repulsion, 
which reduces ^-pT ^^^id makes it coupling-dependent. 
While we employed random-matrix theory, these find- 
ings can be verified for individual systems by varying the 
coupling between their symmetry-related parts. 
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